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C ~ ~) ' The astrophysical S-factor for 16 O + 16 O is investigated within the adiabatic molecular 

0O , picture. It very well explains the available experimental data. The collective radial mass 

C ~ ) ■ causes a pronounced resonant structure in the S-factor excitation function, providing a 

motivation for measuring the 16 O + 16 O fusion cross section at deep sub-barrier energies. 
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1. Introduction 



The fusion cross section a t us at very low energies of reactions involving 12 C and 
16 O is a crucial ingredient to calculate astrophysical reaction rates for different 
stellar burning scenarios in massive stars, in which 16 O + 16 O is the key reaction 
for the later oxygen burning phase. This cross section is usually represented by the 
S-factor (S — <Tf us Ee 27rv , where rj is the Sommerfeld parameter), as it facilitates 
the extrapolation of relatively high-energy fusion data because direct experiments 
at very low energies are very difficult to carry out. Unfortunately, there is a huge 
uncertainty in the S-factor resulted from the extrapolation of different phenomeno- 
logical parametrizations that explain the high-energy data, as shown by Jiang et 
alPfor 16 + 16 0. For 12 C + 12 C, the presence of pronounced molecular resonance 
structures makes it much more uncertain^. These extrapolated values result in reac- 
tions rates that differ by many orders of magnitude^. Therefore, a direct calculation 
of the S-factor at energies of astrophysical interest (< 3 MeV) is essential. 
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We report on an investigational of the fusion reaction 16 O + 16 O within the 
adiabatic molecular picture^, which is realistic at low incident energies. This is be- 
cause the radial motion of the nuclei is expected to be adiabatically slow compared 
to the rearrangement of the two-center mean field of nucleons. In this reaction the 
nuclei are spherical and coupled channels effects are expected to be insignificant 
(the first collective excited state (3~) of 16 is at 6.1 MeV), making its theoretical 
description relatively simple. Furthermore, abundant experimental dataPI exist for 
comparison to the model calculations. 

A basic microscopical model to describe the studied reaction is the two-center 
shell model (TCSM), a great concept introduced (in practice) in heavy-ion physics 
by the Frankfurt schooM We have used a new TCSM^ based on realistic Woods- 
Saxon (WS) potentials. The parameters of the asymptotic WS potentials including 
the spin-orbit term reproduce the experimental single-particle energy levels around 
the Fermi surface of whereas for 32 S the parameters of the global WS poten- 

tial by Soloviev^ are used, its depth being adjusted to reproduce the experimental 
single-particle separation energies^. To describe fusion, the potential parameters 
(including those of the Coulomb potential for protons) have to be interpolated be- 
tween their values for the separated nuclei and the compound nucleus. The param- 
eters can be correlated^ by conserving the volume enclosed by certain equipotential 
surface of the two-center potential for all separations R between the nuclei. 



2. Calculations and discusion 

The adiabatic collective potential energy surface V(R) is obtained with Strutin- 
sky's macroscopic-microscopic method, whilst the radial dependent collective mass 
parameter M(R) is calculated with the cranking mass formulcff^. For simplicity, 
the pairing contribution to the collective potential and radial mass is neglected. 
The rotational moment of inertia of the dinuclear system is defined as the product 
of the cranking mass and the square of the internuclear distance. T he macroscopic 
part of the potential results from the finite-range liquid drop modet^ and the nu- 
clear shapes of the TCSM^. T he m icroscopic shell corrections to the potential are 
calculated with a novel method^. The TCSM is used to calculate the neutron and 
proton energy levels^ Ej&s a function of the separation R between the nuclei along 
with the radial coupling-^ between these levels that appears in the numerator of the 
cranking mass expression, 

„ (8)=! «^ E »f. (1) 

i=i j>A i 1 

Figure []Jl shows the s-wave molecular adiabatic potential (thick solid curve) 
as a function of the internuclear distance, which is normalized with the experi- 
mental Q-value of the reaction (Q = 16.54 MeV). The sequence of nuclear shapes 
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Fig. 1. (a) The s-wave collective potential energy as a function of the separation between the nuclei 
for ls O + le O. The arrow indicates the geometrical contact separation, (b) The radial dependent 
collective mass parameter (in units of nucleon mass mo). See text for further details. 



related to this potential^ is also presented. For comparison we show the Krappe- 
Nix-Sierk (KNS)potentia£El (thin solid curve) and the empirical Broglia-Winther 
(BW) potentiate (dotted curve). Effects of neck between the interacting nuclei, 
before they reach the geometrical contact separation (arrow), are not incorporated 
into the KNS potential. The concept of nuclear shapes is not embedded in the BW 
potential which tends to be similar to the KNS potential. Comparing the KNS 
potential to the molecular adiabatic potential we note that the neck formation sub- 
stantially decreases the potential energy after passing the barrier radius (Rb = 8.4 
fm). It will be shown that the inclusion of neck effects is crucial to successfully 
explain the available S-factor datcP for the studied reaction. 

Figure [T}) shows the radial dependent cranking mass (thick solid curve), whilst 
the asymptotic reduced mass is indicated by the dotted line. Just passing the barrier 
radius, when the neck between the nuclei starts to develop, the cranking mass 
slightly increases compared to the reduced mass and pronounced peaks appear 
inside the geometrical contact separation. For the studied reaction, these peaks are 
mainly caused by the strong change of the single-particle wave functions during 
the rearrangement of the shell structure of the asymptotic nuclei into the shell 
structure of the compound system. In general, the peaks could also be due to 
avoided crossings^ between the adiabatic molecular single-particle states^, which 
can make the denominator of the cranking mass expression |lj very small. It is 
important to stress that the amplitude of these peaks may be reduced by (i) the 
pairing correlation that spreads out the single-particle occupation numbers around 
the Fermi surface, and (ii) the diabatic single-particle motion^ at avoided crossings, 
which can change those populations. For the strongest peak in Fig. [TJa, which is 
located very close to the internal turning point for a wide range of sub-barrier 
energies, only aspect (i) may be relevant as the radial velocity of the nuclei is 
rather small there, suppressing the Landau-Zener transitions. For compact shapes, 
aspect (ii) may lead to intrinsic excitation of the composite system, but this is not 
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Fig. 2. (Color online) The S-factor as a function of the center-of-mass energy for le O + le O. 
The curves are theoretical calculations, whilst the symbols refer to experimental data. The arrow 
indicates the Coulomb barrier of the molecular potential of Fig[T^,. Sec text for further details. 



important here for the calculation of the fusion cross section. Fusion is determined 
by the tunneling probability of the external Coulomb barrier, as explained below. 

Having the adiabatic potential and the adiabatic mass parameter, the radial 
Schrodinger equation is exactly solved with the modified Numerov method and the 
ingoing wave boundary condition imposed inside (about 2 fm) the capture barrier. 
The fusion cross section <jf us is calculated taken into account the identity of the 
interacting nuclei and the parity of the wave function for the relative motion (only 
even partial waves L are included here), i.e., af us = irh 2 /(2fiE) J2l(^^ + !)(! + 
81,2) Pl, where [i is the asymptotic reduced mass, E is the incident energy in the 
total center-of-mass reference frame and Pl is the partial tunneling probability. 

Figure [2] shows the S-factor as a function of the incident energy in the center- 
of-mass reference frame. For a better presentation, the experimental data of each 
setP are binned into AE = 0.5 MeV energy intervals. In this figure the following 
features can be observed: 

(i) the molecular adiabatic potential of Fig. QJ, correctly (thick and thin solid 
curves) explains the measured data, in contrast to either the results obtained 
with the BW potential (dotted curve) or the very recent calculations within the 
Fermionic Molecular Dynamics (FMD) approaclP^] (dashed curve). Since the 
width of the barrier decreases for the molecular adiabatic potential of Fig. [T^,, 
it produces larger fusion cross sections than those arising from the shallower 
KNS and BW potentials. 

(ii) the use of the cranking mass parameter of Fig. QJ> notably affects the low 
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energy S-factor, which is revealed by the comparison between the thick and 
thin solid curves. It starts reducing the S-factor around 7-8 MeV energy region 
and produces a local maximum around 4.5 MeV. At the lowest incident energies 
(below 4 MeV) the S-factor is suppressed by a factor of five compared to that 
arising from a constant reduced mass. The peak in the S-factor is due to an 
increase of the fusion cross section, which is caused by the resonant behavior 
of the collective radial wave function^. 

3. Concluding remarks and outlook 

The adiabatic molecular picture very well explains the available experimental data 
for the S-factor of 16 + 16 0. The collective radial cranking mass causes a relevant 
peak in the S-factor excitation function, although the pairing correlation (neglected 
in the calculation) may somewhat reduce the magnitude of this bump. It is highly 
desirable to have fusion cross sections, measured around 4-5 MeV, to verify the ex- 
istence of this resonant structure in the S-factor. The collective mass surface is very 
important for the reaction dynamics, and its effect on fusion of heavy-ions should 
be investigated systematically. This can be significant for a better understanding 
of reactions forming superheavy elements. Works are in progress for understand- 
ing the molecular resonance structures in the S-factor excitation function of the 
challenging system 12 C + 12 C. 
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